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RESUMEN 

Globulos neutros sumergidos en regiones HII son frecuentemente encontrados 
en regiones de formacion estelar. En este trabajo investigamos la formacion de 
regiones de sombra detras de globulos neutros iluminados por el flujo ionizante 
directo de dos estrellas. Utilizamos una aproximacion analitica, la cual nos permite 
hacer estimaciones de la geometria de las sombras detras de los globulos neutros para 
distintas configuraciones geometricas iniciales, y tambien simulaciones numericas 
tridimensionales. Encontramos un buen acuerdo entre los calculos analiticos y los 
numericos. Una aplicacion de particular relevancia pueden ser los proplyds ubicados 
en el cumulo del Trapecio en Orion, los cuales estan siendo fotoevaporados por la 
radiation de las estrellas O 1 Ori C y 2 Ori A. 

ABSTRACT 

Neutral clumps immersed in HII regions are frequently found in star forma- 
tion regions. We investigate here the formation of tail of neutral gas, which are not 
reached by the direct ionizing flux coming from two massive stars, using both an 
analytical approximation, that allows us to estimate the shadow geometry behind 
the clumps for different initial geometric configurations, and three-dimensional nu- 
merical simulations. We found a good agreement between both approaches to this 
theoretically interesting problem. A particularly important application could be 
the proplyds that are found in the Trapezium cluster in Orion, which are being 
photoevaporated primordially by the O stars 1 Ori C and 2 Ori A. 

Key Words: HYDRODYNAMICS — H II REGIONS — STARS: PRE- 
MAIN SEQUENCE — STARS: FORMATION 



1. INTRODUCTION 

High density neutral clumps immersed in a bath 
of photoionizing and photodissociating radiation can 
be found in several regions in our own Galaxy. To 
quote only a few examples, we can find comet ary 
globules like CG 1 in H II regions (in the Gum 
Nebula), for which the major source of u ltra-violet 
(UV) photons is the Zeta Puppis star fe.g.. lR,eipurthl 
1983); the proplyd s found near the Trapezium clus - 
ter in Orion (e.g., l0 7 Del]lll994 iBallv et allll99fih . 
as well as in other H II regio ns with star for- 
mation (e.g., ISmith et al.l 120031) . The morphol- 
ogy, dynamics and emission line spectrun of the 
Orion proplyds is explained as the photoevaporation 
of circumstellar disks by th e ionizing photon flux 
that comes from 1 Ori C (I Johnstone et al J h 99*1 
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iRichling &; Yorke1l2QQ0l) . Another example of high 
density neutral condensations in an ionizing radia- 
tion field are Thackeray's globules in the IC 2944 
region, which show a fragmented, clumpy structure 
(|R,eipurth et al.ll2QQ3h . We can also find high density 
clumps associated with planetary nebulae (PN). The 
Helix Nebula (NGC 7293, e.g., lO'Dell fc Handronl 
Il996l) displays thousands of cometary shaped neu- 
tral clumps, which are be ing photoevaporated (e.g., 
iLopez-Martm et ai1l200ll) by the radiation field that 
comes from the central star. Attempts have also 
been made to reproduce observational aspects of 
the so-called Fast Low Ionization Emission Regions 
(or FLIERs) observed in a fraction of PN as a by- 
product of the photoevaporation of a neutral clump 
by th e stellar radiation field (e.g., iMellema et al.l 
Il998h . 

Previous analytical studies of high density neu- 
tral clumps exposed to an ionizing photon ra- 
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diation field show that these sys tems have two 
distinct dynamical phases (e.g.: lBertoldi| Il989t 
iBertoldi fc McKeel 1199(1 iMellema et al.lll998l) : an 
initial collapse (or radiation- driven implosion) phase 
and a cometary phase. The collapse phase starts 
when the ionizing photon flux reaches the clump sur- 
face, forming a D-type ionization front (IF). This 
IF is preceeded by a shock front, which propagates 
through the clump (at a velocity determined by the 
local isothermal sound speed) compressing, heating 
and accelerating the clump material. The ionized 
gas expands outwards (with a velocity approximately 
equal to the isothermal sound speed of the ionized 
material) generating a photoevaporated flow 4 . The 
recombination of the gas downstream of the IF (i.e., 
in the photoevaporated flow) partially shields the 
clump surface from the ionizing photon flux that 
comes from the ionizing source, and this is also 
an important effect influencing the evolution of the 
collapse phase. When the shock front propagates 
through the entire clump, in a timescale given by 
to = 2R c /v s (where R c is the initial clump radius and 
v s is the shock speed), the clump starts to accelerate 
as a whole (as a dynamical response to the produc- 
tion of a photoevaporated flow; i.e., the rocket effect) 
and enters the cometary phase. The clump material 
is found to be completely ionized in a timescale that 
depends on both to and t he shielding of the im ping- 
ing ionizing flux (e.g., see IMellema et al.lll998l) . 

One of the first numerical studies of an ionizing 
radiation field interacting with a neutra l clump was 
relate d to t he problem of st a r form ation. iKlein et alJ 
11113) andlftandfordeUD developed a two- 

dimensional code in order to numerically reproduce 
the interaction of an ionizing radiation field from (an 
already formed) OB star association with local inho- 
mogeneities in a molecular cloud. These early stud- 
ies were able to capture the radiation-driven implo- 
sion of the clump material. They have also found 
that the shock induced by the ionization front at 
the clump surface substantially increases the density 
when compared with an analytical, one-dimensional 
evaluation. 

The subsequent evolution of these neutral clumps 
to the cometary regi me was investigated numeri cally 
in two-dimensi o ns by ILefloch fc LazarefJ (|l994l) and 
IMellema et al.l (|l998l) . In particular, the numeri- 
cal simulations of a comet ary globule carried out by 
ILefloch fc Lazarefj (|l994l) show that, after a short 
timescale (~ 10% of the clump lifetime) associated 

4 Actually a neutral clump could be instantly ionized with- 
out evolvi ng through the collapse phase. See a detailed dis- 
cussion m IBertoldil (11989ft . 



with the collapse phase, the clump evolves to a situ- 
ation of quasi-hydrostatic equilib rium, characte ristic 
of the cometary ph ase fsee also | Bertoldilll 9891) . In 
a subsequent paper, ILefloch fc Lazarefj (|l995) show 
that the cometary globule CG7S could be success- 
fully explained as a neutral clump undergoing the 
collapse phase under the influ ence of a nearby group 
of O stars. On the other hand. lMellema et alJ ([19981) 
were able to reproduce the kinematic and emission 
properties of FLIERs seen in PN with a model of 
a clump, located in the outer parts of a PN, being 
photoionized by the central PN star. They have also 
been able to follow the onset of the collapse phase 
and its photoevaporated flow, the effect of the pho- 
toevaporated flow on the clump shaping, and the 
further acceleration of the clump in the cometary 
phase. More recently, three-dimensional numeri- 
cal simulations of non -uniform high density clumps 
([Gonzalez et al .112 0051) sub je ct also to the influence 
of a wind (|Raga et al.l f2005) have been carried out. 
These effects are important for the study of proplyds 
as well as for Thackeray's globules. 

The observed high density neutral structures as- 
sociated with H II regions and PN, project a shadow 
away from t he direction o f the impinging ionizing 
photon flux. ICanto et alJ (|l998l) studied the shape 
and structure of the shadow projected by a spheri- 
cal clump in a photoionized region, also taking into 
account the effect of the diffuse ionizing field. They 
were able to describe (analytically and numerically) 
the transition between a shadow region that is opti- 
cally thin to the diffuse ionizing radiation field to a 
neutral inner core inside the shadow. They also pre- 
dicted that the Ha emission coefficient is substan- 
tially greater in the shadow region when compared 
with the surrounding H II region (and, then, to con- 
clude that the shadow should be brighter than the 
H II region). A sim ilar approach was followed by 
IPavlakis et alJ ()200ll) . who carried out 2-D simula- 
tions. They found that if the diffuse ionizing field is 
10% of the direct field, the evolution of the clump is 
considerably different to the evolution of the clump 
without the diffuse radiation field. Nevertheless, 
none of these papers have attempted to follow the 
whole evolution of the shadow as the clump evolves 
from the collapse through the cometary phases. 

Observationally, some photoionized neutral 
structures ( for example, the Orion proplyds, see 
lO'Delliggci) have relatively short elongations of neu- 
tral material into the shadow region, while others 
have extremely long "tails" extending away from the 
photoionizing source. Examples of long tails are the 
cometary knots in the Helix nebula (see, e.g., the re- 
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cent paper of |Q 'Dell, Henney Ferland||2005|) and 
the dark trails which cut through the outer regions 
of M 42 (jO'DelllteOOnh . In these tail regions, one has 
a combination of the shadowing effect and the pho- 
toionization due to the diffuse radiation (see Canto 
et al. 1998 and Pavlakis et al. 2001) as well as the 
possible presence of a photodissociated wind coming 
from the back side of the neutral clump or disk struc- 
ture (as has been explored for the c ase of the Orion 
proplyds bv lRichling fc YorkelEoonh . 

An interesting effect is that neutral clumps in 
H II regions in some cases are being photoionized 
by the radiation from more than one star. For ex- 
ample, the structures of some of the proplyds ap- 
pear to be affected not only by the radiation from 
1 Ori C, but also by 02 Ori A (in particular the 
197-427, 182-413 and 244-440 proplyds, see O'Dell 
1998 and Henney & O'Dell 1999). A possibly more 
com plex example are th e proplyds in the Carina neb- 
ula (Smii bh et aljEool . which are possibly affected 
by the radiation from several stars in the neighbor- 
hood. This situation was first addressed by Klein, 
Sanford and Whitaker (1983), who performed two- 
dimensional radiation-hydrodynamics calculations of 
the interaction of the radiation field from two mas- 
sive stars with a neutral clump, in order to inves- 
tigate whether or not such a system could trigger 
star formation within OB subgroups. The calcula- 
tion performed by Klein et al. considered two identi- 
cal stars (separated by 1 pc) and, in the middle of the 
straight line connecting both stars, a neutral clump 
(with a radius of 0.6 pc). They found that the clump 
is strongly compressed (by a factor of 170) after its 
implosion and reaches the Jeans mass at some points 
(located in a torus that surrounds the clump initial 
position) in a time scale smaller than that in which 
the clump can be photo-evaporated. Thus, they con- 
clude that stars could be formed in OB associations 
by this " multiple implosion mechanism" . 

In the present paper, we investigate the evolu- 
tion of a high density neutral clump subject to the 
influence of two ionizing photon sources. In some 
sense, this paper represents a generalization of the 
Klein et al. (1983) work, since they only studied the 
axysimmetric case in which the stars and the clump 
are aligned 5 . We first present an analytical descrip- 

5 The assumption made by Klein et al, in which both stars 
are identical and located at the same distance from the clump 
implies that the clump should not be accelerated since it 
should actually be photo-evaporated at the same rate on both 
sides faced to the stars. Also, this situation leads to a mini- 
mum shadow configuration when the distances from the clump 
to the stars are much greater than the clump radius: a situ- 
ation that we will adopt in the present paper. Since one of 
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Fig. 1. Schematic diagram showing i) the relation be- 
tween the two coordinate systems, O and 0\ ii) the 
clump represented by a sphere of radius R c (projected 
onto the yz— plane), iii) the angles a and f3 and iv) the 
distance from the star R*. See the text for a discussion. 

tion of the shape of the shadow behind the clump 
considering both the distance from the sources as 
well as their angular displacement 6 . A set of three- 
dimensional numerical simulations is also presented, 
in which we follow the evolution of the shadows as 
well as the clump radius and neutral mass. 

This paper is planned as follows. In §2, we 
present an analytical solution for the shape of the 
shadows behind a clump illuminated by two sources. 
In §3, we explore the parameters numerically and 
we present the results of three-dimensional numeri- 
cal simulations of this problem. In §4, the discussions 
and conclusions are presented. 

2. SHADOWS BEHIND ILLUMINATED 
CLUMPS 

Figure E] defines two frames of reference O and 
O' . Both frames have their origins at the center of 
a spherical clump of radius R C1 and the star is on 

our goals in the present work is to follow the evolution of the 
shadow behind the clump as it is accelerated by the rocket 
effect, and since that the axi-symmetric case has been already 
treated in the literature, we will not address here the case of 
a aligned system. 

6 In a pioneering study, iDvsonl (jl973ft found that a clump 
subject to an angular distribution of ionizing radiation field 
should respond contracting its radius and increasing its local 
density towards the maxima in the radiation field. Our ap- 
proach differs considerably from this paper since we also fol- 
low the dynamical evolution of the clump and its associated 
shadow. 
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star 1 



star 2 // 



Fig. 2. Schematic diagram showing the relative position, 
angles and distance between the two stars and the clump. 
See the text for a discussion. 

the yz— and y'z'— planes. Frame O has an arbi- 
trary orientation, while frame O' is oriented such 
that its z'-axis is along the line joining the center 
of the clump with the star. The star is located at 
a distance from the clump and has coordinates: 
x = 0, y = —R^s'ma and z = —R*cosa in frame O; 
and x' = 0, y' = and z' = R* in frame O' (see 
Figure CJ. The half opening angle ft of the shadow 
cone is given by: 



sin/? 



R c 
R~*' 



(1) 



Thus the width uo of the shadow cone increases as, 



tan/?(i?* + z) 



Rc 



— (R+ + z >). (2) 



The shadow cone is tangent to the clump at, 

z' c = Rcsinft, y' c = R c cosft, (3) 
thus the shadow region behind the clump is given by, 



(x' 2 + y' 2 ) 1 ' 2 <cu = tan/?(#* + *'), (4) 



and 



z > z' c = — RcSinft, 



(5) 



where the equal sign in (@J corresponds to the bound- 
ary of the shadow. 



^ penumbra 




Fig. 3. Schematic diagram showing the umbra and the 
penumbra regions in the shadow behind a clump illumi- 
nated by two stars. 

The coordinate transformation equations be- 
tween frames O and O' are, 



x 

y' 



-- X, 

y cosa - 
= y sina 



z sma, 
- zcosa. 



(6) 



Substituting (0 in (HJ) and © we find that the 
shadow region in frame O is, 



[x 2 + (ycosa — zsina) 2 ] 1 ^ 2 < 
R 

( R 2 ^2)1/2 + ZC ° Sa + y silia )i ( 7 ) 



and 



y smft 



z cosa > — — - . 

R* 



(8) 



Consider the case of a clump illuminated by two 
stars, one located at a distance R\ at an angle ql\ 
and the other at a distance R2 with an angle 012- 
The yz— plane contains both stars and the center of 
the clump (see Figure EJ). The shadow produced by 
the two stars is divided in two regions. One, the 
umbra where the light from both stars is blocked by 
the clump, and the other, the penumbra where only 
the light from one of the stars is shaded by the clump 
(see Figure 01 . 

In the yz— plane (x = 0) the boundary of the 
shadows are straight lines given by (see eq. 0, 
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ycosai — zsinai = ±tan/3i(i?i + zcosai + ysmai) 

(9) 

for star 1, and, 



ycosa 2 — zsina 2 = ±tdiii/3 2 (R 2 + £cosa2 + ysina 2 ) 

(10) 

for star 2, where, 



tan/3i 



(i?f-i??)l/2 



and tan/?2 



(#1 



i?2)l/2- 
(11) 

Assuming that a 2 > then the maximum ex- 
tent of the umbra Ro (see Figure is given by the 
intersection of lines (0 (with the plus sign) and line 
(fTUI) (with the minus sign). The coordinates of the 
intersection point are, 



xo = 0, yo = Rosinao and = i^ocosao, (12) 
where 

R c 



Ro 



sm7 



and 



a 



(ai +a 2 ) - (f3 2 - Pi) 



(a 2 - ax) - + /? 2 ) 



(13) 



(14) 



(15) 



The geometrical interpretation of ao and 7 is shown 
in Figure [3 

It follows from (|T3|) that Ro —> 00 and 7 — > 0, for 

(a 2 — <^i) (A + /fe)- Therefore for (a 2 — a±) < 
(Pi + /fe) the umbra is not bounded and it extends 
to infinity. 

Let us define a function J(x,y, z,a, R*) such that 
J = 1 if Q) an d @ are satisfied and J = oth- 
erwise. Then, for the umbra J(x, y, z, ai, Ri) = 1 
and J(x, y, z, a 2 , R 2 ) = 1 while for the penumbra ei- 
ther J(x,y,z,ai,Ri) = and J(x,y,z,a 2 ,R 2 ) = 1 
or J(x,y,z,ai,Ri) = 1 and J(x,y,z,a 2 ,R 2 ) = 
0. In the region illuminated by both stars 
J(x, y, z, ai,Ri) = and J(x, y, z, a 2l R 2 ) = 0. 

There is one case of particular interest, that in 
which both R\ and R 2 are much greater than i? c , and 
that we will explore in the next section numerically. 
In the limit, (R c /Ri) and (R c /R 2 ) — > 0, both f3\ and 
f3 2 ^ and thus [see {2J) and (|T5j) ]. 



a 



(ax + a 2 ) 



and 7 



(a 2 - ai) 



(16) 



with Ro given by (|13|) . 

To get an idea of the shape of the umbra in this 
case let us consider a particular orientation of our 
frame of reference. This orientation is such that 



a 2 = 7r — a\ . 



(17) 



Then, from (TL6| . 



ao = I and 7=^-^1, (18) 



and from (p"3|) . 



cosai 



(19) 



The boundary of the shadow projected by each 
star is [see (0], 

x 2 + (ycosai — zsinai) 2 = (20) 
for star 1, and, 



(ycosa 2 — zsina^) 2 = R 2 



(21) 



for star 2. Solving (J20J) and J2U simultaneously and 
using (|T?j) it follows that the condition for the in- 
tersection of the boundaries is z = 0. That is, the 
shadows intersect at the plane z = 0, which is con- 
sistent with the value of ao = tt/2 [see (|T8|) ]. In this 
plane, the shape of the umbra is then, 



y cos ai = R\ 



(22) 



which represents an ellipse centered at the center of 
the clump with semi-axes R c in the x-direction and 
Ro = R c /cosai in the ^/-direction. 

In order to compare the results from these ana- 
lytical solutions with those from the numerical sim- 
ulations, let us anticipate here some results before 
starting to discuss the simulations in detail in the 
next section. Although we have explored different 
angular distributions of these two stars with respect 
to the clump position, let us take here the case in 
which 27 = 45° (see Figure EJ), and a distance from 
the clump to the two stars of d = 3.6 x 10 18 cm. This 
situation corresponds to model M45 discussed in the 
next section (see below). In Figure 21 we show both 
the analytical solution for the shadow geometry (top) 
and the ionization fraction (which should trace the 
shadow behind the clump) obtained through fully 
3-D numerical simulations (bottom) of model M45. 
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t = 30,000 yr 













Fig. 4. Top: A spherical clump and its projected umbra 
and penumbra on the xz- (left) and 2/z-planes (right), 
assuming that the two stars are on the yz-plane. Bottom: 
Maps of the ionization fraction at the beggining of the 
calculation (t = 2 000 years). As in the previous maps, 
a spherical clump creates a shadow on the xz- (left) and 
yz-planes (right). In all cases, the abscissa corresponds 
to the ^-coordinate and is limited by 40 < z < 130, and 
the ordinate corresponds to the x (left) or y (right), both 
limited by 19 < x,y < 109 (in units of 7.8 x 10 15 cm). 

Also, Figure 21 shows the shadows on the xz- (left) 
and yz-plane (right). The stars and the clump are on 
the ?/z-plane. We note that there is a good agreement 
between the predicted geometry for the shadow (top- 
diagrams) compared with the ones obtained from the 
simulations (bottom-diagrams). As the simulations 
progress in time, the clump is photoevaporated and 
assumes a non-spherical shape, which can be seen in 
the bottom-diagrams of Figure [5J Such a departure 
from a spherical shape produces small differences be- 
tween the predicted (see top-diagrams in Figure 0) 
and the simulated shadow geometry. 

3. THE NUMERICAL SIMULATIONS 

3.1. The numerical method and the simulations 

In order to investigate the temporal behavior of 
the shadows of neutral (or partially neutral) material 
behind the clump exposed to two ionizing sources, we 
have carried out a set of three-dimensional numer- 
ical simulations. The simulations were performed 
using the Yguazu-a code. The Yg uazu-a code, a bi- 
nary adaptive grid code (see, e.g.. lRaga et alJEoOOl 
has been extensively used in the literature, 



Fie:. 5. The same as in Figure H but for a more evolved 
temporal stage (t = 30 000 years). In all cases, the ab- 
scissa corresponds to the z-coordinate and is limited by 
80 < z < 170 (as the clump is accelerated towards the 
positive z direction; see §3 below), and the ordinate cor- 
responds to the x (left) or y (right), both limited by 
19 < x,y < 109 (in units of 7.8 x 10 15 cm). 

and has been tested against analy tical solutions and 
laboratory experiments (see, e.g., iRaga et alJl200(j 

1ET 
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200 it IRaga fc Reipurthll2004l) . The Yguazii-a code 
integrates the gas-dynamic equations (employing the 
"flux vector splitting" scheme of Van Leer). The 
code also solves rate equations for neutral/ionized 
hydrogen, and the radiative tra nsfer of the direct 
photo ns at the Lyman limit (see iGonzalez fc Ragal 
l2004> . 

We have computed models assuming a high den- 
sity (n c = 5 x 10 4 ), low temperature clump (T c = 10 
K) immersed in a low density (n e = 100 cm ~ 3 ), 
high temperature environment (T e = 1000 K). All 
the models assume a clump of radius r c = 10 1T 
cm located at 1 parsec from the ionizing sources. 
The ionizing sources were assumed to have a black- 
body spectrum with an effective temperature of T* = 
5 x 10 4 K and an emission rate of ionizing photons of 
S* = 1.5 x 10 49 s _1 . In all models, the computational 
domain is limited to y = 10 18 cm, z = 2 x 10 18 
cm and a 5-level binary adaptative grid with a max- 
imum resolution of 7.81 x 10 15 cm along each axis 
has been used. The centre of the clump is initially 
located &t x = y = z = 5 x 10 17 cm. 

In order to explore the geometrical effect on the 





t=15,000 yr 




Fig. 6. Distribution on the yz— plane of the total (left) and hydrogen neutral densities (right) for model MO (see Table 
QJ at t — 5,000, 10,000, 15,000 and 30,000 years (from top to bottom). The bars on the top give the scale in g cm -3 . 
The ionizing front comes from the left part of the computational domain (the ionizing sources are not depicted here). 



shadows due to the presence of two ionizing sources, 
we have computed models with distinct relative an- 
gles between the clump and the two sources. In Table 
[I] we give the angles 0, which are equivalent to 27 as 
previously defined in Figure [2 (see also equations H5I 
andCnj). In particular = for model M0, = 45° 
for model M45 and = 90° for model M90 (see Table 



[Q . We also note that we have conducted simulations 
with different ratios between the ionizing fluxes from 
the two sources, F = i/S^, namely, F = 1, and 
10 (with = 1.5 x lb 49 s" 1 ) . In the next sec- 
tion, we will discuss in detail the results from these 
numerical experiments. 



CERQUEIRA ET AL. 



10 22 10 20 10 18 10 27 10 24 10 21 10 18 



1 1 III 11 III 1 1 1 1 


1 i i 1 i i 1 i -, | j 


density 

# 

t=5,000 yr 


HI 


• 

t=10,000 yr 








t=15,000 yr 




m 




t=30,000 yr 





Fig. 7. The same as in Figure El but for model M45 (see Table Q). 



3.2. Results 

Figures El and El show the temporal evolution 
of the total density (left) and the neutral hydrogen 
density (HI; right) for models M0, M45 and M90, re- 
spectively (see TableHJ) at t = 5,000, 10,000, 15,000 
and 30,000 years. In these density distribution maps, 
we can identify i) a tail that fills the shadow region 
behind the clump (the umbra in Figure OJ), ii) the 



emergence of a photoevaporated flow (that propa- 
gates towards the ionizing sources), and in) the ac- 
celeration of the clump due to the rocket- effect. 

As discussed in §2, the umbra is expected to have 
an extension Rq, which is in turn controlled by two 
parameters, namely, the clump radius R c and the 
relative angle 7 with respect to the sources (see equa- 
tion El an d Figure EJ). In particular, the simulations 
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Fig. 8. The same as in Figure El but for model M90 (see Table Q). 



show that the tail which appears in the shadow re- 
gion behind the clumps has a cylindrical shape for 
model M0 {Rq — > oo; Fig. EJ, while in models M45 
and M90 the umbra depicts a conical shape (see Fig- 
ures Hand [HJ). The conical umbra, however, is some- 
what smaller in model M90 when compared with the 
one of model M45. We also note that the the cone 
height Rq is actually a function of time {Rq = R{i)), 



since the photoevaporation of the clump material re- 
duces the clump radius (see the discussion below). 
In all cases, both the umbra and the penumbra are 
completely symmetric with respect to the z-axis. 

The photoevaporated flow can be seen in both 
the total density and HI density maps in Figures 
El El and [HI It expands almost radially from the 
clump position with a velocity < 20 km s _1 , pro- 
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Fig. 9. The same as in Figure El but for model M45b (see Tabled). 



ducing density perturbations in the environment (a 
shock wave). This photoevaporated flow accelerates 
outwards, towards the ionizing photon sources. As a 
consequence, the clump itself is accelerated in the op- 
posite direction. This effect is clearly seen in Figures 
H3 and [SI where the clump is continuously pushed 
towards larges distances from the ionizing sources. 
However, we also note that the mass flux in the pho- 



toevaporated flow is smaller in model M90 than in 
model M45, and, as a consequence, the acceleration 
of the clump is also lower (see the discussion below) . 

In Figure EH we show the evolution of the total 
density (left) and the neutral hydrogen density (HI; 
right) for model M45b (see Table [TJ at t = 5, 000, 
10,000, 15,000 and 30,000 years. We note that in 
this particular model, in which the ionizing sources 
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TABLE 1 
THE SIMULATED MODELS 



Model 


on 


Si is' 1 ) 


^2 is' 1 ) 


MO 





1.50 x 10 49 


1.50 x 10 49 


M45 


45 


1.50 x 10 49 


1.50 x 10 49 


M45b 


45 


1.50 x 10 49 


1.50 x 10 48 


M90 


90 


1.50 x 10 49 


1.50 x 10 49 



have different ionizing photon rates (S^i/S^ = 10), 
both the photoevaporated flow and the shadows be- 
hind the clump, particularly the penumbra, are not 
symmetric with respect to the x— axis. In particular, 
the photoevaporated flow is stronger in the direc- 
tion of the most powerful ionizing source, Si (lo- 
cated towards the top-left direction of the diagram; 
see Figure EJ. This causes a more pronounced ab- 
lation of the clump material that is facing the more 
powerful ionizing source. Also, only the penumbra 
of the shadow produced by the ionizing field from Si 
is seen. 

In FigureslTUk andlTHb we show, respectively, the 
position Zcm of the center of mass of the clump and 
its velocity as a function of time (in years), for mod- 
els M0 (full line), M45 (dotted line), M45b (dashed 
line) and M90 (dot-dashed line) 7 . As mentioned be- 
fore, the clump in model M0 is rapidly accelerated 
(when compared with the clump in the other models) 
and attains large distances from the ionizing sources 
at shorter times. This is due to the fact that the 
mass loss rate associated with the photoevaporated 
flow, M c , is higher in this case. This is illustrated 
in Figure ITUb . where the clump mass (normalized 
by the initial clump mass, M c (t = 0) = M Cy o = 
po(4:7r/3)Rl = 4.55 x 10 32 g) is shown for all mod- 
els, as a function of time. Note that the clump mass 
loss in the M0 model is higher (when compared with 
the other models). 

In Figure HUH we also show the spherical clump 
radius R c (normalized by the initial clump radius, 
R c (t = 0) = R Ci o = 10 1T cm), as a function of 
time, for all of the models listed in Table All 
of the curves show an almost monotonic behaviour. 
Note that the clump radius in model M0 is always 
smaller than the values obtained from the other mod- 
els (which is consistent with the higher M c of model 

7 We note that the coordinates Xqm and Yqm are al- 
most constant, and equal to the initial value Xcm{^ = 0) = 
Ycm it = 0) = 5 x 10 17 cm, for all models with the excep- 
tion of model M45b, where Ycm changes with time. Hence, 
the match in the curves from model M45b (dashed line) and 
M90 (dot-dashed line) in this diagram does not mean that the 
clump has the same velocities. 




^ 1.0 - 



0.5 - " 

2x10 4 4x10 4 6x10' 

Time (years) 



Fig. 10. a) Clump centre of mass (CM) position (mea- 
sured by Zcm, in units of 10 17 cm) as a function of 
time for models M0 (full line), M45 (dotted line), M45b 
(dashed line) and M90 (dot-dashed line) (see Table HJ. 
We note that the variations in Xcm and Ycm are almost 
null for all models with the exception of M45b, see the 
discussion in the text, b) Velocity of the clump centre 
of mass (Vcm'i in units of km s _1 ) as a function of time 
(the curves are labeled in the top panel), c) Clump mass 
(normalized by the initial clump mass, M c ,o = 4.55 x 10 32 
g) as a function of time (the curves are labeled in the top 
panel), d) Clump radius (normalized by the initial clump 
radius, R c ,o = 10 17 cm) as a function of time (the curves 
are labeled in the top panel), e) Shadow parameter Ro 
(see eauation HSJ) . normalized by the initial clump radius, 
as a function of time (we note that Ro — > oo for model 
M0; not shown here). 

M0). The temporal variation of the clump radius is 
responsible for the temporal variation of the height 
of the conical shadow i?o, as can be seen in figures 
□ M and M In Fig. HQt, the height of the conical 
shadow Ro (normalized by the initial clump radius, 



12 



CERQUEIRA ET AL. 



25 




2x10 4 4x10 4 6x10 4 

Fig. 11. Top: Clump centre of mass (CM) position (mea- 
sured by Zcm, in units of 10 17 cm) as a function of time 
for simulated model MO (full line; the same as in Figure 
HU|) . Also depicted in this figure is the res ult from the an- 
alytica l approximation (diamonds) from Mel lema et alJ 
( 1998), using as input parameters the sound speed Co and 
Ci taken from the numerical simulation (see the text). 
Bottom: Clump mass (normalized by the initial clump 
mass, M C; o = 4.55 x 10 32 g) as a function of time taken 
from the MO model (full line; the same as in Figure EOt 
. As in the previous diagram, the diamonds represents 
the analytical solution. 

Rc,o', Ro —> oo for model MO, not shown in Figure 
ITUbl is calculated 8 using the clump radius given in 
Fig. fTUll and equation JT3J). Note that, after a strong 
variation at the beginning (t < 10 4 years), the conical 
height tends to values « Ro(t = 0)/2 in all models. 

We have also compared the results from our sim- 
ulations with results obt ained from the analy tical 
approximations given by Me llema et al.l {l998). In 
that paper, the authors studied both the collapse 
(or implosion) and the cometary phases of a clump 
being exposed to an incident plane-parallel ionizing 
front (the case of model MO here), and the mass 
and the position of the clump as a function of time 
have been derived for both phases. Using values ob- 
tained from our simulations for model MO (following 
the notation in Mellema et al.: sound speed at the 
base of the photoevaporated flow, q ~ 8.64 km s _1 ; 
sound speed of the shocked gas inside the clump, 
Co — 2.12 km s _1 ), we found that the collapse phase 
takes to = 13 234 years. Interestingly, almost all the 
curves in Figure EH seems to change their behavior 
at this time, and this seems to be due to the end of 

8 These curves show the qualitatively behavior of Ro with 
time, since we are calculating it at each time step, without 
taking into consideration the adjustment of the shadow to the 
steady- state. 



this phase and the beginning of the cometary phase. 
At time to, the clump is predicted to have a mass 
M c ~ 0.797 x M c>0 . Figure ITOh shows that model 
MO has a mass at this time that is ~ 10% smaller 
than this value, in good agreement with the analytic 
theory. Finally, the behavior of the solid curves in 
Figure ITUk.c has a lso been compared wi th the solu- 
tions presented bv lMellema et all (|l998l) and a very 
good agreement is found between the values from the 
simulations and those predicted by the analytic the- 
ory during all of the evolution of the clump in model 
MO. 

In particular, in Figure ITT1 we plot again the re- 
sult from model MO (full lines) and its comparison 
with the result from the a nalytical approximation 
from iMellema et al.l (1 19981 ) (diamonds). Although 
the clump CM position is well recovered in the sim- 
ulation (top panel), it seems that the mass of the 
clump is under-estimated in the numerical simula- 
tion when compared with the analytical solution. 
However, even in this case, both curves (full line and 
diamonds, on the bottom panel of Figure ITTj) depicts 
the same behavior. We have also simulated this same 
model MO, but with a grid-resolution of 512x128x128 
(these numbers corresponding to a uniform grid at 
the highest resolution of our adaptative grid), or a 
grid spacing of 3.9 x 10 15 cm (i. e., an improvement 
of a factor 2 with respect to the simulations that we 
have presented up to now). As in the previous, low 
resolution case, the CM position is well reproduced 
by the numerical simulation (not shown here), and 
the same under-estimation of the clump mass is also 
obtained. 

4. CONCLUSIONS 

In this paper, we have explored analytically what 
is the shape of the shadow behind a clump exposed 
to two ionizing photon sources. For the case in which 
the clump radius is much smaller than the distances 
to the sources, an analytical solution for the shapes 
of the umbra and the penumbra are found. We have 
also carried out 3D numerical simulations of this sce- 
nario, and we present four models with different rela- 
tive positions between the two stars and the clump as 
well as with different combinations of ionizing pho- 
ton rates for the sources. 

These models show that the long neutral tails 
produced in the interaction of a neutral clump with 
a single photoionizing source (model MO, see Table 
1 and Figure 4) disappear for the case of two an- 
gularly separated photoionizing sources producing 
comparable ionizing photon fluxes at the position of 
the clump (models M45 and M90, see Figures 5 and 
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6). This shortening of the neutral tails is again ob- 
tained in a model in which the ionizing photon fluxes 
from the two sources differ by a factor of 10 (model 
M45b, see Table 1 and Figure 7). In our numerical 
simulations, we also obtain the time-evolution of the 
cometary clump as it is accelerated away from the 
photoionizing sources by the rocket effect. 

The present models are meant as an illustra- 
tion of the effects of the interaction of a neutral 
clump with two photoionizing sources, which might 
be applicable to neutral structures in H II regions 
with several massive stars (see, e.g.. ICVDelll llQQfit 
ISmith et al.l l2003h . However, in our simulations 
we have not included the diffuse ionizing radia- 
tion (which might have an important effect on the 
shadow region, in particular in the case of ionization 
bounded nebulae) and the effect of the dissociative, 
direct and diffuse radiation (whic h can affect both 
the h ead and the tail regions, see iRichling h Yorkel 
l2000l) . In order to model the structure of specific 
clumps, these effects should also be included. 
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